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Abstract 

We show that the infinite percolating cluster (with density Poo) of central- force networks is 
composed of: a fractal stress-bearing backbone (Pb) and; rigid but unstressed "dangling ends" 
which occupy a finite volume-fraction of the lattice (Pd). Near the rigidity threshold p*, there is 
then a first-order transition in P^o = Pd + Pbi while Pb, is second-order with exponent (3' A new 
mean field theory shows (3'^^ — 1/2, while simulations of triangular lattices give (3[ = 0.255 ±0.03. 
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Figure 1: The geometry of the infinite rigid cluster, a) An infinite rigid cluster at the rigidity 
threshold of a site-diluted triangular lattice (p* = 0.69707 for this configuration). Backbone 
(stressed) bonds are black, dangling ends (rigid but unstressed bonds) are blue and cutting bonds 
(if they are removed the backbone collapses) are red. b) Poo (long dashes); Pb (solid lines) and; Pd 
(short dashes) for site- (plusses) and bond- (circles) diluted triangular lattices. The calculations 
were done at the percolation point for each lattice size and are averaged over roughly 10^/L^ 
configurations. 

The connectivity-percolation geometry has become a paradigm in the study of disordered 
systems[l], especially close to the percolation critical point. However, in many problems involv- 
ing mechanical properties, such as the mechanical properties of granular media[2], glasses[3] and 
gels[4], the connectivity-percolation geometry does not apply. In systems such as these, in which 
central forces are of primary importance, the infinite cluster must be multiply connected in order to 
transmit stress and hence to support any mechanical property. Here we present a comprehensive 
analysis of the geometry of the infinite cluster in this "rigidity percolation" problem. Early work 
on the stress-bearing paths in central-force systems relied on direct solution to the force equations 
for model systems [5-8], such as lattices of Hookian springs. More recently, a more accurate method 
which finds the backbone[9] and infinite-cluster[10] geometry directly has been developed. There 
has also been progress in the development of mean-field theories, both continuum[ll] and using 
Cayley trees, and we report on the latter here. Using a combination of these new techniques a 
novel picture of rigidity percolation has emerged. In particular, we show that: 

• The infinite cluster is composed of a stress-carrying backbone and rigid (but unstressed) 
"dangling ends" (see Fig. la). The backbone is fractal and this is the reason that the elastic 
constants undergo a second- order transition on approach to the rigidity transition [5-9]. 

• The "dangling ends" occupy a finite volume fraction of the lattice and for this reason there 
is a first- order transition in Poo at 

• Due to the fractal backbone, there is a length which diverges with a non-trivial exponent, 
V, at the rigidity threshold. This exponent is different than the connectivity-percolation 
correlation-length exponent. 

• The number of "red" or cutting bonds scales with exponent 1/v aX the rigidity threshold. 

Consider a triangular lattice composed of nodes connected by Hooke springs. Consider attach- 
ing rigid "bus-bars" on two sides of the lattice, and then removing 1 — p of the sites (or bonds) 
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Figure 2: a) One branch ofsiZ = 4 Cayley tree. The lowest-level branches of the tree are attached 
to a rigid border, b) P^o for site-diluted trees with g = 2 and, starting with the rightmost curve, 
z = 4,6,9,13. 

from the lattice. Upon application of an external stress (e.g. a tensile stress) to the lattice, we 

find results such as that presented in Fig. la. In this figure, the black bonds carry stress, while 
the "blue" bonds are rigid but do not carry stress. The red bonds carry stress and are critical in 
the sense that if they are removed, the lattice can no longer support the applied stress (they are 
"cutting" bonds) [12]. The configuration in this figure is exactly at the percolation point. This 
figure was found using a new "graph-theory" method[13] which enables us to find the rigidity- 
percolation geometry exactly and which has been described elsewhere [9, 10]. We plot the volume 
fraction of backbone, dangling-end and infinite-cluster bonds as a function of sample size in Fig. 
lb. It is seen in this figure, that although the volume fraction of "dangling" bonds is rather small 
(roughly 0.1), it is constant, indicating that the infinite-cluster probability undergoes a first-order 
jump at the rigidity threshold. The stressed backbone however shows a non-trivial scaling, so that 
Ub ~ I/^^, with Db = 1.78 ±0.02 as found previously [9, 10]. Note that a fit to the infinite-cluster 
data of Fig. lb suggests that the infinite cluster is fractal, as claimed by Thorpe and Jacobs[10]. 
However, by separating the dangling ends from the backbone bonds as shown in Fig. lb, the 
weakly first order character of Pqo is evident. Thus in the thermodynamic limit, we predict, 

Poo=a + b{p-p,f, (1) 

where for the triangular lattice, the backbone exponent (3f = 0.255 ± 0.03. Here we have used 
the result, 2 — fi' /u = Db, with u = 1.16±0.03 from our previous calculations [9, 10]. For the 
site-diluted triangular lattice the first-order jump in Pqo, a*^ = 0.086 ± 0.005, while for the case 
of bond dilution = 0.11 ± 0.02 (see Fig. lb). Our mean-field theory uses exact "constraint 
counting" (see below) on trees. We assign each node of a lattice g degrees of freedom and co- 
ordination z. On the triangular lattice treated in the previous paragraph, each node has two 
degrees of freedom and 2; = 6 for the pure lattice. If the nodes of a triangular lattice were extended 
objects or bodies instead of being pointlike "joints", then each node would have an additional 
rotational degree of freedom. Cases of physical interest are then: 

5 = 1 for connectivity percolation (2) 

g = d for a joint (3) 

g = d + '^^^ ^ ^'^ for a body, (4) 

where d is the dimensionality of the embeding space. Knowledge of g and z is enough to write down 
the simplest constraint-counting mean- field theory [3]. Consider bond dilution with probability p 
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that a bond is present. Present bonds restrict the motion of the nodes, and hence they impose 
"constraints" on the degrees-of-freedom of the lattice. We assume that all of the constraints are 
"independent" (a "dependent" bond can be removed without causing any reduction in the size of 
the rigid clusters), then the number of degrees of freedom that remain unconstrained (the "floppy" 
modes per site /) is approximated by, 

f = g- pz/2 for p<p^^^. (5) 

For p > p^, / = in this approximation. This mean- field theory has been useful in the study of 
glasses [3], but has yielded little information about Poo, Pb and Pp. More recently, a continuum 
mean field theory has been developed [11], which focuses on Pqo- That theory indicates that the 
infinite-cluster probability undergoes a first-order transition at the rigidity threshold, but the 
connection with the key parameters g and z are unclear. In addition, in that work a pathological 
lattice model(a square lattice with random diagonals) was suggested as a paragidgm for the rigidity 
transition. We show later why that model is anomalous. 

We now present our Cayley-tree model for rigidity percolation, which provides a complete mean- 
field model for arbitrary g and z. The trees have co-ordination number z but, as usual in tree 
models, the key results come from consideration of one branch of a tree (see Fig. 2a). Consider 
site-diluted trees which are grown from a rigid boundary at infinity. Building inward from this 
boundary, we keep track of the number of degrees of freedom of a node with respect to the 
boundary. Rigidity can only be transmitted to higher levels of the tree if there are enough rigid 
bonds present to offset the g of degrees of freedom of a newly-added node. For connectivity 
percolation only one bond is needed. If a node is added to a 51 = 2 tree, two bonds are needed 
to offset the two degrees of freedom of the added node. In general, if the nodes of the tree have 
g degrees of freedom, rigidity is transmitted to the next level of the tree provided the node is 
occupied, and provided at least g of the lower-level nodes are rigid. This gives the recurrence 
relation, 

k=g 

Where is the probability that a site which is I levels from the rigid border (which is level 0) 

is rigid with respect to the border. For g = I this reproduces the familiar model for connectivity 
percolation[14], while for 51 > 1 it is equivalent to so-called "bootstrap" percolation[15] (Note that 
the equivalence of the transmission of rigidity and bootstrap percolation does not apply to regular 
lattices). If we take the thermodynamic limit (very large /), Eq. (6) iterates to a steady-state 
solution, which we call T^a- Finally we find Poo from 



iT^)'{l-T^r-\ (7) 



The results for Too for several values of g and z are presented in Fig. 2b. It is seen from this 
Figure that for g = 2, the rigidity transition is first order, we show elsewhere that this conculsion 
holds for all g > 1[16]. An exact solution is straightforward for the case z = 4, g = 2, where (from 
Eq. (6)) 

T^=piTl + 3Tl{l-T^)), (8) 
which has the trivial solution Too = and the non-trivial solutions 

= (<,) 

In order to ensure that T^o = 1, when p = I, take the positive solution in Eq. (9). The square 
root in Eq. (9) becomes imaginary at = 8/9. As p — > p*, Eq. (9) yields, 

Too =3/4 + c(p-p.)i/2 (10) 
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Figure 3: The volume fraction of red bonds Pr at the rigidity threshold for: The triangular 
lattice with site dilution (inverted triangles) and bond dilution(triangies); The "random diagonal" 
model (see text) with q = (squares-this is the anomalous model of Fig. 4), g = 0.1 (circles) and 
q = 0.4(diamonds). 

The key feature of this equation is that Too is first order and there is a square-root singularity 
superimposed upon the first-order transition. We find a similar behavior for all g and z, and so 
we write in general, 

P^=a + b{p-p,y/\ (11) 

where a, b and are dependent on g and z. Notice that this is of the same form as Eq. (1) for 
the triangular lattice, and that the square-root behavior is due to the backbone. We thus argue 
that in general the form Eq. (1) is correct, with the mean-field backbone-exponent = 1/2. 
Full details of the Cay ley-tree result will be published elsewhere [16]. 

Due to the fractal backbone, there is a diverging length, and the scaling properties are controlled 
by this diverging length. As shown in previous calculations this "correlation" length diverges with 
exponent ly = 1.16±0.03 [9,10] at the rigidity threshold in two dimensions. Coniglio's relation then 
applies to this problem [12], so that we expect Ur, the number of red or "cutting" bonds to diverge 
as L^^^ at the rigidity threshold. We find that to be the case for both the site- and bond- diluted 
triangular lattices (see Fig. 3), from which we find P,- ^ Ur/L"^ ^ L^~^, with x = 0.86±0.02 = l/i^ 
confirming our previous result [9]. Finally, we analyse and extend a simple model for gels, which 
has been developed by Obukov[ll] to provide an illustration of a first-order rigidity transition 
(see Fig. 4 - the boundary conditions do not change the conclusions). The random diagonals 
are present with probability pd, and all that is required to make the lattice rigid is one diagonal 
present in every row of the square lattice. The probability that an infinite-cluster exists is then 
P_l_ = (1 — (1 — Pd)^)^ , from which we find pd* ^ InL/L. However, the resulting infinite cluster 
contains the whole lattice, the rigid backbone is extensive, and so is the number of cutting bonds 
(see Fig. 3 for a calculation of the cutting bonds). This is inconsistent with the result Eq. (1) and 
also with the presence of a diverging length at the rigidity transition. However, the model of Fig. 
4 is atypical, as can be seen by considering a generalized model in which we randomly add the 
diagonals (with probability pd) to a square lattice whose bonds have been diluted with probability 
q Obukov's model is g = 0), while ii q = 1 — pd this model is equivalent to the bond-diluted 
triangular lattice. We find that even for a small amount of dilution, e.g. q = 0.10, the rigidity 
transition returns to the behavior characteristic of the isotropic triangular case (see Fig. 3). We 
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Figure 4: The square lattice with random diagonals and no removed horizontal or vertical bonds, 
so g' = (see text). 

find that for sufficiently large lattice sizes, the universal behavior Eq. (1) holds for any finite 
q < 0.5, and that the "fully-first-order" transition (i.e. a first-order backbone) only occurs in the 
special case of a perfect square lattice (or other "marginal" regular lattices) with randomly added 
diagonals. We have found similar pathological cases on Cayley trees. 
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